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00 ' We work out the thermodynamic equations for the decays and scat- 

■ terings of heavy Majorana neutrinos including the constraints from 

Q_i' unitarity. The Boltzmann equations depend on the CP asymmetry 

parameter which contains both, a self-energy and a vertex correction. 
^ ■ At thermal equilibrium there is no net lepton asymmetry due to the 

CPT theorem and the unitarity constraint. We show explicitly that 



Abstract 



deviations from thermal equilibrium create the lepton asymmetry. 
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1 Introduction 



Over the past few years it has been shown that a lepton asymmetry can be 
generated by the mixing of heavy Majorana neutrinos. Majorana neutrinos 
have the remarkable property of being mixed states of particles and antipar- 
ticles. By definition they induce AL = 2 and consequently A{B — L) = —2 
transitions. In addition they could have couplings to scalar particles, which 
allow them to decay into Higgs particles and leptons, i.e. — > (f)H,(j)i'^ 
In these models CP— violation is introduced through complex couplings and 
appears in the self-energy ^ or vertex corrections |]I|, ^ |, ^. We shall 
classify the effects using the terminology of mesons We shall call di- 
rect or e'— type effects, the CP— violation which arises from the interference 
of tree diagrams with vertex corrections. Similarly we call indirect or (5— type 
effects, those which arise through the self-energies. The self-energies play the 
role of the box diagram in the system and define the physical states. A 
consequence of the observation in 0, ^ is that the physical Majorana states 
\E' are not CP eigenstates and that the CP— asymmetries produce observable 
effects. Finally, the mass splitting between the states also plays a role and 
when the mass difference is of the order of the width there is a resonance 
enhancement 0. 

The processes under discussion must satisfy unitarity constraints so that 
the sum of the probabilities for all transitions to and from a state i should 
sum to one and yields 0, |^: 



For our case this means that when we consider the scatterings icj)^ i'^cf) 
equation (|l|) turns into 



(1) 



j j 



(2) 
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so that the probabihties for the direct and inverse process should be equaL 
For the early universe a summation over all possible physical states with 
Boltzmann factors weighting the inital states is required. Boltzmann fac- 
tors are introduced according to the thermal properties of the universe. At 
thermal equilibrium, for example, the weighting factors are all equal, which 
implies that the lepton asymmetry averages to zero 0. This issue and its 
implications are explicitly proved in P] and further discussed in [|T0|] . 

The generation of a lepton or baryon asymmetry is a combination of the 
Sakharov effects [jll| (i) baryon or lepton violation, (ii) C and CP— violation 
and (iii) the thermodynamic properties of the ensemble of particles in the 
early universe. For this reason several early papers studied the non-equilibrium 
equations of the ensemble |T2|, |r^. These thermal considerations together 
with several subtleties of the phenomenon are mentioned in [Q]. In this ar- 
ticle we present an explicit calculation of the modification of the unitarity 
constraints as implied by the thermodynamics of the system and point out 
the meaning of the various terms. In fact we shall show that the s-channel 
physical poles play a special role and are responsible for the development of 
asymmetries. 

To begin with we consider massive unstable particles which interact and 
decay. As the particles interact some of them occur in intermediate states. 
When the lifetimes of the particles are long relative to the interaction times, 
the intermediate states interact many times with the thermal bath. For this 
reason they should not be considered as virtual states, but as ensembles 
of particles with their own thermodynamic distributions. This phenomenon 
modifies the unitarity sum, i.e. the various terms in equation @ are weighted 
by different factors, and a lepton asymmetry is created. 

We estimate the collision of the particles in relation to their lifetimes. 
The width of the particles with energy E and mass Mn is given by 




for E = T 



(3) 
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Their frequent interactions are like e~ + t N + b with a Higgs particle 
exchanged in the t-channel. The cross-section at temperatures T = E ^ Mm 
is 

We now calculate the interaction of one particle N with a density of quarks. 
At a temperature T ^ and assuming Boltzmann statistics the particle 
density is given by 



n=AT3 . (5) 



Therefore the interaction rate for relativistic velocities is 



|2 



n-<^-v= T . (6) 

Thus there are many more interactions per decay width 



with the last factor being of the order of 10® or larger. Since many scatterings 
are taking place during the lifetime of the particles, each state has a differ- 
ent history and the particle states are incoherent. At this epoch when we 
consider many scatterings, decays and recombinations of particles, there is 
no asymmetry, because the unitarity of the S-matrix guarantees the equality 
of direct and inverse processes. An asymmetry is generated when specific 
channels begin to decouple. 

As usual we work on an extension of the standard model where we in- 
clude one heavy right-handed Majorana field Ni per generation of light-lepton 
{i = 1,2,3). The new fields are singlets with respect to the standard model 



|l|. Extensions of this model have been constructed in gauge theories |14, 15 



supersymmetric theories |T^, and theories with electroweak singlet neu- 
trinos |T^ . 

We consider the contribution from the interference of the diagrams in 
Figure 1, which appears in the asymmetry of the reaction £^0 i(f)^. In 
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Figure 1: Interference of the tree with the one-loop self-energy diagram 



the figure, we have indicated the momenta explicitly and have taken the 
absorptive part of the self-energy. In addition we consider many generations 
indicated by the indices a and f3 for the external particles, the index 7 for 
the loop-momentum and i,j and k for the heavy particle propagators. The 
contribution to the asymmetry is 

2Re[MLM) - ^^,{kp + pp kp) _ Mf){s - Ml){s - M]) 

- ^^'iky I , . M.M,Im(/^,,/^.,/^^/^^,/^,,/.;.) 

- ^(fcp + PP kp) (8) 

with AIl is the loop diagram and M. is the Born diagram amphtude (see 
Figure 1). Carrying out the same calculation for the process l(f^ — > I'^cf) we 
obtain the following expression for the asymmetry: 

2lie{Mr^M)- ^^,{kp+pp kp) _ Mf){s - Mi){s - Mf) " 

In this expression we kept the propagator field for all intermediate states 
but omitted their widths. The two expressions are the same except for the 
factor of the coupling constants. A careful comparison of the two expressions 
shows that they become identical when we sum over the intermediate states 
and over all external particles. This can be seen by rearranging of the in- 
dices. It is a consequence of unitarity mentioned above. The same property 
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persists when we keep all the masses of the propagators or when we replace 
the propagators with 5— functions in the narrow-width approximation. The 
two expressions 2Re(A12-^) ^-iid 2Re(Al|,A^) are different when we select 
specific components. This is in fact what happens in the development of the 
universe, as the temperature reaches the mass of a heavy neutrino specific 
components are multiplied by special Boltzmann factors, which single them 
out as demonstrated below. 



2 The Boltzmann Equations 

In this section we explicitly discuss the thermodynamic development of the 
system. In this way we explain how deviations from thermal equilibrium 
modify the sum of unitary contributions and render them different from zero. 
The lepton asymmetry is a sum of several terms: 

First, there are the decays of the heavy Majorana neutrinos: 

— > which creates an excess of leptons, and 
"^i — > £'^(f) which reduces the amount of leptons 

and the corresponding recombination terms. In addition there are two scat- 
tering processes: 



£'^4' (.(fy which creates an excess of leptons, and 
£(j)^ £'^(f) which reduces the amount of leptons. 

The result of all these terms is the following differential equation for the 
lepton asymmetry density = — n^c : 

riL + SHriL = hA\M{^i^i(j)^)\^ -\M{^i^r(j))\^]Al2 (10) 



+ 



-hU \M^<I'^ ^ ^i)l' + UU\M{r<p ^ * 



or 
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where Afl is the four-particle phase space factor depending on the momenta 
Pi to Pi, fe and /^t are the Boltzmann distributions for leptons and scalar 
particles, respectively, fi^ is the distribution for the heavy Majorana states. 
Since the \E''^s are produced copiously through the multiple scatterings as 
discussed in the previous section they have their own distribution. For the 
decays and recombinations we introduce the three-particle phase space A^2- 
The first two lines already include the decays of the real intermediate states 
(RIS) and their recombinations. For this reason they are subtracted from 
the scattering amplitudes. 

The interference of the tree with the one loop graph gives the CP— vio- 
lating factor: 



5)\Mo\ 



with |A^oP being the tree level amplitude. Both e' and 6 are CP violating 
parameters with e' is produced by the vertices (direct CP— violation) and 
6 from the self energies. Combining the terms and neglecting those ones of 
0{[6' + 6]^ , [fi{e' + 6)] /T), the first two lines of equation (jT^) give: 

n^, {e' + 6) {T^,) + n^f {e' + 6) {T^,) - 



T 



j^w /r \ 



(12) 



where the phase space was combined with the amplitude |A^oP to produce 
the thermally-averaged decay width (F^^). 

Next we consider the last two lines of equation (p!0|), where we expand 
the particle distributions in powers of (yu/T) : 



2 A^^ 



+fiUt\MRis^^^^ 

9 A 34 rEQ 



(13) 
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+ {1 + ^)\Mrjs 



m'-ii-^)\MRjs{t<p 



+0 



J^2 



The leading terms proportional to unity involve the difference of the com- 
plete amplitudes and vanish by virtue of unitarity. This is the only 
place where unitarity is effective in bringing a complete cancellation. The 
remaining terms arise from interactions which change the number of par- 
ticles (chemical potentials) and through deviations from equilibrium. The 
former is the origin of terms proportional to the chemical potentials which 
are multiplied with CP— conserving amplitudes. The leading term from the 
real intermediate states is of special interest, because it is twice as big as 
the second term in equation ([T^), but has the opposite sign. Its net effect is 



to change the sign of the second term in equation ([T2|). After some algebra, 
described in the appendix, equation (O) leads to: 



(14) 



-2 



{a'{i(f)^ ^0) ■ v) + {a'{t(f) £0^) ■ v) 



We can substitute the chemical potential by the density for the lepton asym- 
metry using the equation 

''^ . (15) 



— = 2^ + ^ 



T 



2^3 



Combining the various terms we arrive at the final Boltzmann equation: 



is' + S) 



EQ 



1 

2 



n 



*?(r^.) (16) 



v) 



This is coupled with the differential equation for the density of the Majorana 
neutrinos, which reads: 



+ SHuy^^ = — (r^^)(?2^j 



n 



1 
2 



^(^)<r(^' + 5)(r^,) . (17) 



EQ 



"1 . 
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The integration is now straight forward. We can integrate equation(|T7p, 
then introduce the answer in equation (|TB]) whose numerical integration gives 
the lepton asymmetry. The terms proportional to Ul in equation (|T^) are 
multiplied by negative coefficients and after integration, if they were alone, 
they would wash out any lepton asymmetry exponentially. 



3 Numerical Integrations 

Following the usual method to evaluate the coupled Boltzmann equations 



[jT3| we introduce the dimensionless quantity 

, (18) 

to change the integration variable from the time t to which has a tempera- 
ture dependence. Furthermore it is convenient to look at the particle density 
per comoving volume, which is: 

F = ^-^ , (19) 

with s being the entropy density and g^, the number of effectively massless 
degrees of freedom (particles with mass m <^T). Now the differential equa- 
tions ([iBD and (13) read: 

dYL _ (rvl>i(2 = 1))^ f^^, , ^^^^^ ^^EQ. 



^ = -%^^fif--.{i'*.-n?(i+o(.'.^))} (20) 

with 

(T^^iz)) Ki{z) { z z<^\ 
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(r*,(^ = i)) K^{z) 11 ^»i 
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where are the modified Bessel functions |]^, and 

{\g^Y^^{V^^{z))+2{a'v)n,] 



^3/2e-^ + O.lz"^ z>l 



Here we approximate the cross-section terms as in |T^. For a detailed calcu- 
lation of the cross-section terms see 0. 

The overall factor 

--'^ 

is an important measure for the efficiency to create a lepton asymmetry. In 
the case K <^1, which corresponds to (r^^(z = 1)) <^ H{z = 1), the decay 
rate is much smaller than the expansion rate of the universe and the particles 
come out of equilibrium and create a lepton asymmetry. In the other case, 
where K ^ 1, the particle decay fast and recombine and are in thermal 
equilibrium in comparison to the expansion rate of the universe. In this case 
the lepton asymmetry will go to zero. 

Figure 2 shows this dependence of the lepton asymmetry, which we obtain 
by integrating the Boltzmann equations (|20|) numerically. The shape of the 



curves is in direct correspondence to |]13l, but the asymptotic value of the 



curves is now bigger, because of the CP— parameter 6 coming out of the 
interference of the tree with the self energy diagram: 

Yl{z ^oo) = , for fsT < 1 (22) 

9* 

If K increases the asymptotic value of the lepton asymmetry decreases and 
finally for K ^ 1 goes to zero. Similar numerical results were obtained in 
[1T9| where the energy scale for the generation of the lepton asymmetry was 
studied. 
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K=0.01 




0.01 0.1 1 10 100 1000 

z = Mi/T 

Figure 2: The development of the lepton asymmetry Yl/ [(e' + 5)/ g^}\ with 
the expansion of the universe given hy x — Mi/T 



4 Conclusions 

We have shown explicitly how a lepton asymmetry is created through the 
decay of Majorana neutrinos. The asymmetry is proportional to the sum of 
the self energy and the vertex contribution. We emphasize that at thermal 
equilibrium, the unitarity of the amplitudes implies the vanishing of the 
lepton asymmetry, as it is demonstrated in section 1. On the other hand, 
the thermal development of the Majorana densities is the reason for the 
unitarity cancellation to become ineffective and allow the generation of the 
lepton asymmetry (see section 2). Finally, we solve the Boltzmann equations 
numerically and create a final asymmetry, which is proportional to the sum 
of the CP— parameters coming from the vertex contribution e' and from the 
self energy S. The latter one has a resonant behaviour depending on the mass 
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difference of the Majorana particles 
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A The Mixed Majorana States 

The problem of mixing of Majorana neutrinos has been treated at several 
places. One formalism is to discuss the development of the various com- 
ponents in the wavefunction, where the absorptive part has a definite effect 
P, 0]. Another method is to construct the inverse propagator and their renor- 



malization |T^, |T^ . In both cases the absorptive part is a physical observable, 
which can not be removed by renormalization |T^. We follow here the ap- 
proach of references 0, where to 0{h?^^ the CP— effects were computed 
exactly for all relative values of the masses 0. 

We assume a mass hierarchy and consider an epoch in the development of 
the universe, where the Majorana neutrinos are decoupled from each other. 
When the temperature of the universe becomes of the order of the lightest 
mass Ml, there is only this particle left, which can only decay bringing the 
universe out-of-equilibrium. The physical heavy Majorana Neutrino \E'i is a 
mixed state of the interaction states: 

l^i) = i= {\Nr,) + \{Nr{)-) + «2|A^ffi) + «i|(A^R2)^)} (23) 



It has a definite mass and equation (^3]) gives the couplings of the various 
components. The transition amplitudes of the mixed state receive contribu- 
tions from the vertex corrections through the e' and from the self energies 
through the mixing of the states: 
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I-Mc 



:i + e') + 



(24) 



— k:7^ — + 



(25) 



with I^VloP = J2a \hai\'^Mf. These are the relevant amphtudes for the Boltz- 
mann equation (|l3), which are given now by: 



+fi' 



I^Qll 



Re(a2h*^ha2) 



(26) 



Re(h;^h^2(a2 - «i))' 



+2 Alt {UU [\M{t^ ^ £0t)|2 _ \Mnis{^'^ - 

Next we substitute the phase space densities and develop them in the pa- 
rameter /i/T : 



EQn/T 



EQ 



1^1 



EQ-^/T 



EQ 



(27) 



As usual we assume that the particles are in kinetic equilibrium. Finally, we 
use the definition of the thermally-averaged decay width (F) which is: 

g/{2nY J d'pfip)T{p) g/{27^f J d?pf{p)V{p) 



(r) 



(28) 



gl{27:f!d^pf{p) 

where g are the degrees of freedom and n is the particle density. With all 
this information we arrive at equation(|T2|) and equation (13). In the next 
section we study the scattering terms in the Boltzmann equation. 
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B The Scattering Amplitudes 



In this section we study the scattering terms in the Boltzmann equation: 



+ 



fx 



\M' 



(29) 



with 



\Mf = \M\'-\MRis\' (30) 
For the real intermediate states one can use the narrow-width approximation: 



TT 



\MRis{t 



TT 



Using the following integral |H]: 



5{s - Ml)\M{t<p ^ ^i)\''\M{^i 



(31) 



r 1 

J dU,dU,f^^6{s - Ml) = (^2vr^^^(r^,) (32) 
and the definition of the thermally-averaged cross section 

gA/i27Tf9B/i27Tf I d^PAl d^PBfiPA)fiPBMAB ^ CD)v 



{a-v) 



QaI {27^^93 /{27Tf J d^pA I d^PBf{pA)f{PB) 



(33) 



gA/{2nfgB/{2nY I d'pAl d^PBf{pA)f{PBMAB CD)v 



UaUb 



as well as the approximations ~ n^c n^t ^ n-y, expression (|T^ can 

be derived from equation (p9|). 
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